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Chirality is one of the key features governing the electronic properties of single- and bilayer graphene: the 
basics of this concept and its consequences on transport are presented in this review. By breaking the inversion 
symmetry, a band gap can be opened in the band structures of both systems at the K-point. This leads to inter¬ 
esting consequences for the pseudospin and, therefore, for the chirality. These consequences can be accessed by 
investigating the evolution of the Berry phase in such systems. Experimental observations of Fabry-Perot inter¬ 
ference in a dual-gated bilayer graphene device are finally presented and are used to illustrate the role played by 
the band gap on the evolution of the pseudospin. The presented results can be attributed to the breaking of the 
chirality in the energy range close to the gap. 


I. INTRODUCTION 

Experimentally isolated in 2004 [1], single-layer graphene 
(SLG) consists of a layer of carbon atoms, arranged in a hon¬ 
eycomb pattern. Its unit cell is defined by two carbon atoms, 
usually referred to as A and B, forming the two-atom basis of 
a Bravais lattice. 

As predicted in 1947 [2], it was experimentally demon¬ 
strated in 2005 that charge carriers in graphene behave like 
massless Dirac Eermions [3, 4]. They can be described by a 
two-component wavefunction yi obeying: 

- ihvY^{G • V)i/r = Exfi, (1) 

where vp is the Eermi velocity and a = (cTx, (Ty) is a vector 
of two Pauli matrices. Here, the analogy with quantum elec¬ 
trodynamics can be made by realizing that the two sublattices 
A and B play the role of spin-up and spin-down and that G is 
not the spin but the pseudospin operator. The direction of mo¬ 
tion is coupled to the pseudospin orientation, as one can see 
from Eq. (1), a property denoted as chirality. The chirality of 
charge carriers has important consequences for transport. It is 
responsible for a Berry phase of ;r [3, 4] and the suppression 
of backward scattering [5, 6]. These two effects are the basis 
of the Klein paradox [5, 7]. 

Bilayer graphene (BLG) also exhibits chiral charge carri¬ 
ers. However, instead of following the Klein physics, BLG ex¬ 
hibits anti-Klein properties, due to a Berry phase of In [5, 8]. 
The concept of chirality in both single- and bilayer graphene 
and its presence in interference experiments is the focus of 
this review. 

In Section II, we introduce the concept of chirality in SLG 
and BLG, and lay an emphasis on illustrating this concept, to¬ 
gether with the concept of Berry’s phase. We then focus on 
several interference experiments, where signatures of the chi¬ 
rality were successfully observed. In Section III, we consider 
SLG and BLG systems in which the inversion symmetry has 
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been lifted and a band gap has been opened in the band struc¬ 
ture. We show that this lifting results in strong out-of-plane 
perturbations of the pseudospin in ^-space in the energy range 
close to the gap, but that the pseudospin orientation is further 
restored to its original state at higher energies. Einally, we 
focus on dual-gated BLG in Section IV and present an exper¬ 
iment allowing for the observation of the consequences of the 
opening of a band gap on the chirality, probed in a Eabry-Perot 
interferometer geometry. 


II. PSEUDOSPIN AND CHIRALITY 

Charge carriers in graphene are chiral. This means that 
there is a handedness of their states because their pseudospin 
is locked to the direction of motion, giving rise to interesting 
tunneling properties that we review in this section. 


A. Pseudospin motion in SLG and BLG 


In graphene, low-energy charge carriers live around two 
inequivalent points in momentum space, K and K', called 
“valleys”. In the vicinity of the K-point, single- and bilayer 
graphene can be described by the Hamiltonians [9, 10] 


tisLG = Vp 



^BLG 


^ / 0 (7r'^)2\ 

2m* 0 ) ’ 


( 2 ) 


respectively, where 7t = p^^ipy is the momentum operator 
and m* is the effective mass. The above effective Hamilto¬ 
nians are related to the tight-binding models through vp = 
^y^ajlh with 7o and a the intralayer nearest neighbor hop¬ 
ping energy and distance, respectively, and nf = y\l2vp with 
7i the interlayer nearest neighbor hopping energy. 

The Hamiltonians (2) act on the spinors (v^a,V^b)^ and 
(V^Bi 5 WtsiY ^ respectively, where A and B are the two inequiv¬ 
alent carbon sites of a single-layer graphene and the indices 1 
and 2 refer to the top and bottom layers of bilayer graphene; 
see Pigs, la and lb. This gives rise to the band structure shown 
in Pig. 2, left panels. Note however that the Hamiltonian (2) 
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is a simplified two-band version, well-suited to describe the 
system at low energy. The full description would give rise to 
4 bands: two lower bands touching at the ^-point (seen in the 
figure) and two split bands, split away from the lower ones 
by the energy 0.39 eV, not represented in the figure. In the 
following, we consider that we are in the low energy range, 
where only the lower bands are filled, as shown in Fig. 2. 

The associated eigenstates are: 


Wto) - ^ (^^1 


|v4g) - ^ (^^1 


— 


( 3 ) 


where the ± signs are related to the two eigenenergies E = 
EfiVY^k and (j) = arg (kx + iky) characterizes the direction of the 
wave vector A: = (kx, ky) measured from the ^-point. To better 
visualize the details of the motion of the pseudospin which 
is closely related to the Berry phase, we next introduce the 
polarization vector P as a convenient quantity. 

In quantum mechanics, the polarization vector of a spin- 
1/2 quantum state \\f/) = cos |,sin |)^ is given by the 
expectation values of the Pauli matrices 


^{V/|cTx|V/)^ 

P= ( iwlOylw) 1 = 


' sin 0 COS 
sin 6 sin ( 
cosG 


( 4 ) 


where 6 G [0,;r] is the polar angle and (j) G [0,2;r] is the az¬ 
imuthal angle [11]. This vector of length 1 describes the spin 
orientation of | if/) and can conveniently be represented on the 
Bloch sphere [12]. Represented on such a Bloch sphere, the 
quantum state | iff ), which is a superposition of the two basis 
states |—) = (0,1)^ and |+) = (1,0)^, together with its polar¬ 
ization vector are sketched in Fig. 1(c). 

In the case of SLG and BLG, the polarization vector can 
be calculated by replacing | iff) in Eq. (4) with the eigenstates 
from Eq. (3), leading to 


'coscj)^ 
PsLG= I sin0 


( cos20' 

^BLG= I sin20 


( 5 ) 


which describe, instead of the real spin, the pseudospin direc¬ 
tion of Iv/slg) and Iv^blg). respectively, and are particularly 
convenient for visualizing the motion of the pseudospin when 
the momentum rotates. The polarization vectors of Eq. (5) 
have no z-component, due to the absence of diagonal terms in 
the Hamiltonians in Eq. (2). Thus both Pslg and Prlg depict 
a pseudospin restricted to the equatorial plane of the Bloch 
sphere (0 = ;r/2), as sketched in Fig. 2. 

As indicated in Eq. (5), the pseudospin for the case of SLG 
rotates as fast as its wave vector [see Fig. 2(a)], while for BLG 
it winds twice as fast [see Fig. 2(b)]. In both cases, the process 
is only momentum-dependent but energy-independent. This is 
highlighted by the different constant energy cuts in each dis¬ 
persion [middle panels of Figs. 2(a) and 2(b)], for the case of 
the conduction band. For the valence band, on the other hand, 
the pseudospin is inverted with respect to the K point. This 
is referred to as a chirality -fl for electrons (the pseudospin 
is always parallel to the wave vector) and —1 for holes (the 
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FIG. 1. Lattice orientation of (a) SLG and (b) BLG in Bernal stack¬ 
ing. (c) A quantum state | i/r) and its associated polarization vector 
(red needle) represented on the Bloch sphere, (d) When the polariza¬ 
tion vector undergoes an adiabatic evolution at constant polar angle 
0, it encloses a solid angle Q, which describes a portion of the sphere 
here highlighted by a red shading. 


pseudospin is always anti-parallel to the wave vector) [13]. 
Since the two valleys exhibit opposite handedness, the reverse 
situation happens around the K' point, with chirality —1 for 
the electrons and +1 for the holes. 

An important quantity closely related to the pseudospin mo¬ 
tion is the Berry phase [14]. The Berry phase, also known as 
a “geometrical phase” or “Pancharatnam phase”, is a phase 
acquired by a system during an adiabatic cyclic evolution. 
Upon adiabatic evolution along a closed path, the polariza¬ 
tion vector defines a portion of the Bloch sphere and subtends 
a solid angle Q. An exemplary motion is shown in Fig. 1(d), 
where the polarization vector evolves along a circle of con¬ 
stant latitude at polar angle 6. In this case, the solid angle is 
given by fl = Jq Jq^ sin0d9d(j) = 2;r(l —cosO). The Berry 
phase d>Berry, the loop integral of the Berry connection, 
was found to be related to this quantity by [14-16] 


^Berry = ^ 


( 6 ) 


The Berry phase is half the solid angle subtended by the pseu¬ 
dospin during its motion. In case of an evolution in the equa¬ 
torial plane (6 = 7r/2), as in the SLG case, the Berry phase (6) 
is equal to n. In BLG, however, since the pseudospin rotates 
twice as fast as the momentum, the Berry phase is equal to 271 
[see right panels of Fig. 2(a)-(b)]. 
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Focusing on the left interface of the barrier, we can see that 
a right moving electron (belonging to the right branch of the 
conduction band in the region L) will find a perfect match on 
the other side of the interface (C region), since a right mov¬ 
ing hole exhibits the same pseudospin orientation, as high¬ 
lighted in Fig. 3(a) with the encircled pseudospins. Such a 
situation implies a high transmission probability through the 
barrier. The exact solution to this problem has been investi¬ 
gated by Katsnelson et al. [5] in 2006 and described in de¬ 
tail in Ref. [19]. For an electron wave of energy E incident 
on a barrier of height Vq > 0, with translational invariance 
along the transverse direction y, and holes populating the C 
region (Vq > E + /zvF|^y|), the resulting refiection coefficient 
was found to be: 

r = x 

sin cp — ss' sin 'd- (8) 

cos {(p 'd) cos {(p — t^)] —lisinqjcd' 


FIG. 2. (a) Left: Band structure of SLG at the ^-point with its pseu¬ 
dospin orientation at constant energy in the conduction band. Mid¬ 
dle: The pseudospin orientation in space is momentum-dependent 
but energy-independent. Right: The pseudospin motion on the Bloch 
sphere when the momentum is rotated by Iti. The enclosed portion 
of sphere is highlighted with a red shading, (b) Same content as (a) 
but for BLG. 


Hence, the Berry phase (in units of 7l) is equal to the num¬ 
ber of rotations of the pseudospin when the wave vector com¬ 
pletes one rotation around the A: = 0 point. This is why the 
Berry phase is often called the pseudospin winding number 

[17] . This integer number represents the degree of chirality 

[18] . Interestingly, this description is valid as well for multi¬ 
layer graphene. A /-layer graphene system can be described 
by the Dirac-like Hamiltonian: 


Hj = gj 


0 ) 


g\ = Vf 

g2 = vllY\ 

' 83 = 4/ri 


(7) 


The corresponding Berry phase is accordingly /;r [18]. 


B. Transmission through a potential barrier 

The pseudospin, as described above, leads to interesting 
properties when considering a charge carrier incident on a po¬ 
tential barrier. Consider the situation depicted in Fig. 3(a): 
a particle, moving from left (L region) to right (R region) in 
SLG, is incident at energy £’f > 0 on a potential barrier (C 
region) of width d. The barrier of height Vq > £’f is further 
assumed to result in a np'n junction (the prime referring to the 
C region), shifting the charge neutrality point as schematically 
displayed in Fig. 3(a). If the electron propagates ballistically 
in the regions away from the edges of the barrier potential, 
such a geometry resembles a Fabry-Perot interferometer. 


where 5' = sgn(£’) and s' = sgn{E — Vq) are the sign functions, 
(p and t3- are the incident and refractive angles, and q is the 
transmitted wave vector [5, 19]. This enables the calculation 
of the transmission probability, given by / = 1 — |rp. From 
this expression, one finds that normal incidence always leads 
to full transparency of the barrier (r\(p=T}=o = 0), known as the 
Klein tunneling in SLG. 

Experimentally, such a graphene np’n junction can be real- 



FIG. 3. Schematic of a potential barrier engineered in (a) SLG and 
(b) BLG. The charge neutrality points are effectively shifted such that 
the Fermi level lies at the conduction band in the outer regions (L and 
R) and the valence band in the central one (C). The orientations of 
the pseudospins are indicated with arrows (red for electrons and blue 
for holes). 
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ing as sketched in Fig. 3(b). 

For non-zero incidence angles, some ‘magic angles’ appear, 
where the transmission increases sharply to one. The resulting 
conductance, which is proportional to the transmission inte¬ 
grated over all the incident angles, is therefore much smaller 
than in the SLG case. This is a property of practical interest 
in the sense that electrostatic barriers in BLG are then highly 
efficient to confine carriers. For wider barriers, the number of 
resonances increases quickly, as seen already by comparing 
the two cases shown in Fig. 4. 

One should mention at this stage that these results of T{(p), 
based on the Dirac equation, can be reproduced by the tight- 
binding-model-based Green’s function approach, which can 
easily handle arbitrarily shaped barriers and also allows for 
more general band structures [20] . This will be used in later 
sections to implement the band gap in BLG. 


C. Signatures of chirality in interference experiments 


(a) (b) 

FIG. 4. Calculated transmission probability T as a function of the 
incident angle (j) across a np'n junction with fixed carrier density 
^out = 5 X 10^^ cm“^ outside the barrier and fZin = —3 x 10^^cm“^ 
inside the barrier for SLG (red lines) and gapless BLG (blue lines) 
with barrier thicknesses of (a) d = 100 nm and (b) d = 400 nm. 
The given densities Wout and rii^ correspond to the incident en¬ 
ergy Fp ~ SOmeV and barrier height Vq ^ 280meV for SLG, and 
Ef ^ 16meV, Vb ~ lOOmeV for gapless BLG. 

ized by electrically controlling the carrier density inside and 
outside the barrier region. Considering densities ^out = 5 x 
10^^ cm“^ outside the barrier region and ^in = — 3 x 10^^ cm“^ 
within, we show two examples of T{(p) with barrier thick¬ 
nesses of J = 100nm,400nm in Fig. 4 (red curves), based on 
Eq. (8). These density values correspond to an incident energy 
of £’f ~ SOmeV and a barrier height of Vb ~ 280meV. In these 
plots, we can see that the transmission probability at normal 
incidence is one, remains very high for small incident angles, 
and then exhibits additional full-transmission resonances at fi¬ 
nite angles (sometimes) called ‘magic angles’. 

The same problem in the case of BLG was also treated 
in [5, 19] in the absence of the gap. The angle dependence 
of T can be similarly obtained by solving the transmission 
problem. For a direct comparison to SLG, we consider the 
same and ^out that lead to Ep « 16meV,Vb ~ 100meV 
and show again two examples with barrier thicknesses d = 
lOOnm, 400nm in Fig. 4 (blue curves). Unlike for SLG, where 
massless Dirac fermions are always perfectly transmitted at 
normal incidence, a perfect refiection^ at normal incidence is 
observed. This phenomenon is known as anti-Klein tunneling, 
and can be understood in terms of lack of pseudospin match- 


^ This is true, however, only when d ^ Xp where Xp is the Fermi wave¬ 
length. In the case of d <Xp, evanescent modes can lead to finite trans¬ 
mission even at normal incidence. 


As explained above, the transmission across a barrier in 
graphene exhibits a unique behavior due to the linear disper¬ 
sion and the chirality as above-defined, in contrast to non- 
chiral particles [5]. However, a direct measurement of the 
angle-dependent transmission T{(p) was so far not accessible 
[21, 22], since transport experiments usually measure the total 
conductance G involving all contributing angles. This is why 
a true hallmark of Klein scattering has been sought. In 2008, 
Shytov et al. [23] realized that a signature of Klein tunneling 
in ballistic pn'p Fabry-Perot SLG cavities should appear in 
the magnetic field dependence of the interference pattern. 

In the absence of magnetic field, the directly transmit¬ 
ted and twice refiected waves^ [such as the case sketched in 
Fig. 5(a)] interfere with each other with a kinetic phase differ¬ 
ence d>wKB (WKB standing for Wentzel-Kramers-Brillouin), 
which is proportional to the cavity size d and the longitudi¬ 
nal component of the wave vector k^. Under the infiuence 
of an external magnetic field B, the resulting Lorentz force 
bends the semiclassical electron trajectories, some of which 
form closed loops [such as the case sketched in Fig. 6(a)], 
wrapping finite areas and therefore an Aharanov-Bohm phase 
<I>ab, which can be shown (when B is weak) to ho oc B^/kp, 
where kp is the Fermi wave vector. At the same time, B 
also reduces kx as a consequence of the bending and the con¬ 
servation of kp = {k^ X- kyY^^, where ky is proportional to 
B because of the mechanical momentum p ^ p ^ eA with 
the chosen gauge^ A = (0,v5,0) for B = V xA = (0,0,B). 
Thus the increase of B leads to an increase of <I>ab but de¬ 
crease of d>wKB, while the increase of kp leads to a decrease 
of <I>ab and increase of d>wKB- When measuring the con¬ 
ductance as a function of magnetic field and density k^), 


^ Waves of multiple reflections also contribute to the Fabry-Perot interfer¬ 
ence but are only of minor importance, even in ideally ballistic transport; 
see, for example, [24]. 

^ Note that when applying the periodic boundary condition along the trans¬ 
verse dimension (y), this gauge is practically the only choice in order to 
keep the system y-independent [25]. 
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FIG. 5. (a) Trajectory of a wave in real space at zero magnetic field 
incident from region L toward region R across a central cavity C, 
assuming an imperfect transmission. In C, the refiected (dashed ar¬ 
row) and transmitted (solid arrow) waves are slightly offset for clar¬ 
ity. (b) The corresponding trajectory in momentum space with ky = Q 
(slightly offset for clarity) encloses the origin and is equivalent to a 
rotation of the momentum by 2%, resulting in the contribution of the 
Berry phase. 


FIG. 6. For SLG, the situation depicted in Fig. 5 does not exist due 
to the perfect transmission of Klein tunneling. Instead, a finite mag¬ 
netic field is needed to obtain enough backscattering in the region C. 
An example of the resulting trajectory is sketched in (a), (b) The cor¬ 
responding trajectory in momentum space forms a closed loop that 
encloses the origin. At this magnetic field value, the Berry phase of 
71 is picked up, leading to a phase shift in the interference pattern. 


the two competing phase terms result in the Fabry-Perot in¬ 
terference fringes (each stripe corresponding to a constant 
^WKB + ^ab) dispersing with increasing B towards higher 
density in a parabolic pattern [6, 23, 25-32]. 

The cyclotron bending contributes to the phase in a differ¬ 
ent way. At 5 = 0, the trajectory propagating at normal in¬ 
cidence will be fully transmitted. Thus there exists no such 
trajectory as the one depicted in Fig. 5(a) which, in momen¬ 
tum space, would directly enclose the origin [see Fig. 5(b)]. 
In the presence of a finite (and typically weak) B, however, 
closed loops such as that shown in Fig. 6(a) can form because 
of finite backscattering at oblique incidence, and the corre¬ 
sponding momentum-space trajectories will enclose the origin 
as sketched in Fig. 6(b). Due to the existence of a singularity 
at the ^-point (vertex [15, 33]), the ;r-Berry phase is picked 
up. In the measurements, this will therefore lead to a shift 
of the oscillation pattern by half a period. Hence this ;r-shift 
constitutes a direct evidence for Klein tunneling, through the 
analysis of backscattering. 

The observation of the phase jump has been made for SLG 
cavities of various qualities and widths [6, 27-31]. In 2009, 
Young and Kim presented their pioneering experiment and 
reported on the observation of conductance oscillations in a 
graphene bipolar heterojunction [6]. The cavity was electro¬ 
statically induced by the use of a narrow top gate (~ 20 nm), 
allowing for the observation of ballistic interference, as the 
mean free path was estimated to be larger than 100 nm. The 
signature of Klein tunneling and of its perfect transmission 
was further demonstrated by investigating the magnetic field 
dependence of the conductance oscillations, where the ;r-shift 
of the oscillations was seen. The specificity of the measured 
shift of the oscillations agreed with the theoretical prediction 
[23] and was later reproduced qualitatively [26] and quantita¬ 
tively [25] by transport calculations. 

In 2013, a big step was made in the quality of the produced 
devices. Rickhaus et al. [28] and Grushina et al. [29] reported 
on the use of local bottom gates in suspended devices to engi¬ 
neer smooth, electrostatically defined pn junctions with stun¬ 
ning qualities, enabling for the observation of Fabry-Perot in¬ 


terference on length scales larger than 1 pm [34] . Once again, 
the characteristic Berry phase shift was present in the mag¬ 
netic field data. In Ref. [28], the finesse of the cavity and 
the resulting visibility of the oscillation pattern were studied 
in great detail and analyzed in view of the smoothness of the 
electrostatic landscape. 

In the case of BLG, however, the Berry phase of 271 could 
not be identified in a similar way in an interferometer geom¬ 
etry, since a phase of 271 is equivalent to a phase of 0. Un¬ 
til today, the only measurable contribution of this phase was 
observed in quantum Hall measurements, as experimentally 
demonstrated for the first time in Ref. [8]. Moreover, an elec¬ 
trostatically controlled BLG pn’p junction normally requires 
a dual gated cavity where the inversion symmetry is broken 
such that the BLG is no longer gapless. The consequence of 
the band gap on transport and its relation to the pseudospin 
are to be discussed in the following section. 


III. BAND GAP AND PEUDOSPIN 

The pseudospin in SLG and BLG is restricted to the x-y 
plane as a consequence of the lattice inversion symmetry. If 
this symmetry is broken, the situation is different. 


A. Inversion symmetry breaking 

By applying a different potential to the two sublattices, the 
inversion symmetry can be lifted. In SLG, this can be done 
by aligning the flake on a hexagonal boron nitride substrate 
[35-38]. The closely similar lattice structure of graphene and 
hexagonal boron nitride (h-BN) results in locally aligning the 
atomic site A of graphene with a boron atom and its B site 
with a nitrogen atom (or vice-versa). The two carbon sites 
therefore experience different potentials, leading to the bro¬ 
ken inversion symmetry. In BLG, setting the two layers to 
different potentials, using for example external gates [39-42], 
also breaks the inversion symmetry. The experimental real- 
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ization of the BLG symmetry breaking will be treated in more 
detail in Section IV. 

With the inversion asymmetry, the resulting Hamiltonians 
are generalized from Eq. (2) to the form 


^SLG = 



^BLG = 


/ u 

\ 2m* 



(9) 


where u is the asymmetry parameter [39]. 

One observes in Eq. (9) that adding an asymmetry adds di¬ 
agonal terms to the Hamiltonians, i.e. a CT^-component, which 
gives rise to a z-component of the polarization vector itself. 
In this situation, the pseudospin is not bound to the equato¬ 
rial plane any more. This is shown schematically in Eig. 7(a) 
for the case of the conduction band of SLG at the valley K. 
We see that for small momenta (i.e. small energies, close to 
the band gap), the pseudospin points completely out of plane 
(black arrows). On the Bloch sphere, this means that the pseu¬ 
dospin stays aligned with the z-axis and therefore encloses 
no area, giving rise to a zero Berry phase. This is shown in 
Eig. 7(b) with the black needle. The chirality is broken. Mov¬ 
ing towards higher energies, the pseudospin tends to asymp- 



\4Ja) 


\^a) \^a) \^a) 



\^ b ) |j//e) |j//e) 

(c) (d) (e) 

FIG. 7. (a) Pseudospin orientation along constant energy contours 
of the conduction band of gapped SLG (dashed circles), around the 
^-point. The process is now energy-dependent. The pseudospin mo¬ 
tion on the Bloch sphere is shown for four cases, (b) The black needle 
corresponds to the black pseudospin polarization vectors close to the 
band gap in (a): no area is enclosed and the Berry phase is there¬ 
fore 0. (c)-(e) When going away from the gap, the z-component is 
decreased and the enclosed area grows [colored needles correspond 
to the pseudospin vectors sketched in (a)], until chirality is restored 
at higher energies and the pseudospin returns to the equatorial plane 
[red arrows in (a) and red needles in (e)], recovering to the Berry 
phase of K. 





FIG. 8. Berry phase ^pBerry dependence (neglecting trigonal warping) 
on momentum k for an asymmetry u = 50meV applied between top 
and bottom layer of BLG. As shown here, close to the gap (/: = 0), the 
Berry phase is zero. Going towards higher momenta, and therefore 
towards higher energies, the Berry phase recovers its pristine value 
of 2k. 


totically recover its in-plane motion [red arrows in Fig. 7(a) 
and red needles in Fig. 7(e)] and therefore a ;r-Berry phase. 
The chirality is slowly restored [38]. In-between the two ex¬ 
treme cases, the pseudospin is partially z-polarized, leading to 
a Berry phase varying between 0 and K, as shown with the 
blue and green arrows in Fig. 7(a) and the blue and green nee¬ 
dles in Fig. 7(c)-(d). 

Exactly the same happens with BLG when a band gap is 
opened. Close to the edges of the valence or conduction band, 
the pseudospin is fully z-polarized, leading to a Berry phase 
of 0 [43]. When the energy increases, the chirality is slowly 
restored, until the pseudospin returns to the vy-plane and the 
Berry phase is set back to 2k. This is illustrated in Fig. 8, 
where the Berry phase is calculated as a function of momen¬ 
tum (which can be translated into an energy), and will be fur¬ 
ther explained and investigated in Section IV. 


B. Gapped BLG pn junction 


We next focus on the effect of the broken chirality induced 
by opening the band gap on the transmission probability of an 
electron wave incident on a BLG pn junction. As illustrated 
earlier, the perfect reflection across a bipolar potential step in 
BLG is expected [5] for pristine BLG in the absence of the 
band gap. However, only limited literature exists addressing 
the question of the effect of a gap on the tunneling properties 
(see, for example. Refs. [17, 44]). Using the same method as 
Ref. [20], we apply the Green’s function method based on a 
tight-binding model associated with the periodic Bloch phase 
to illustrate the influence of the asymmetry u on the angle- 
resolved transmission T{(p). This asymmetry results in open¬ 
ing a band gap in the band structure of BLG [39-42]. 

We consider transport through an ideally sharp np June- 
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FIG. 9. (a) Top: An ideal BLG w;? junction with step height Vr — Vl = 
lOOmeV and Ef —Vl = 50meV with the asymmetry varying from 
u = 0 (black curves) to u= lOOmeV (red curves). Bottom: Angle- 
resolved transmission T{(p) with the colors corresponding to those 
of the band structures shown in the top panel, (b) Top: Transmission 
at normal incidence T{(p = 0) from (a) as a function of u. Bottom: 
The corresponding Berry phase. The data points at w = 70meV rep¬ 
resented by ★ correspond to 7(0 = 0) ^ 1 (top) and ^ 

(bottom), i.e., a revival of the Klein tunneling in BLG. 


dispersion, but also to tune its chirality, enabling a possible 
switching between BLG-like and SLG-like transport charac¬ 
teristics. In the following, we present an experiment illustrat¬ 
ing how the chirality in a gapped BLG system can be broken. 


IV. INTERFERENCE EXPERIMENT WITH GAPPED BLG 


As mentioned earlier in Sec. IIC, there is no report in the 
literature showing explicit signatures of the 2;r-Berry phase 
in BLG interferometers, as there is in SLG. However, as de¬ 
scribed in the previous section, once a band gap is opened, the 
chirality of BLG is changed and a finite transmission at nor¬ 
mal incidence, associated to a finite Berry phase, occurs. In 
the following, we investigate the chiral properties of BLG via 
measurements performed on a dual-gated BLG Fabry-Perot 
(FP) interferometer. Parts of the data presented below were 
published in Ref. [24]. 



tion^ with a step height of lOOmeV at fixed Fermi energy 
Ef = 50meV as shown in the top panel of Fig. 9(a), where 
the asymmetry parameter u is varied from 0 (black curves) to 
lOOmeV (red curves). The resulting T{(p) curves for various 
values of u in this range are shown in the bottom panel of 
Fig. 9(a), where one observes that the transmission function 
changes drastically. Starting from u = 0 [black curve with 
r((p = 0) = 0], we observe that increasing u breaks the anti- 
Klein tunneling, and T{(p = 0) increases slowly, until reach¬ 
ing perfect transmission at normal incidence (curve marked by 
★). The T{(p = 0) values and the corresponding Berry phase 
as a function of u are respectively shown in the top and bottom 
panel of Fig. 9(b), where T{(p = 0) = 1 and <I>Berry = ^ match¬ 
ing each other at = 70meV can be seen. This implies that, 
under certain circumstances, a revival of the Klein tunneling 
in BLG is possible, by manipulating the gap-controlled Berry 
phase. 

By further increasing u, the Berry phase ^Berry continues 
to increase toward In while the normal incident transmis¬ 
sion r((p = 0) decreases rapidly toward zero, until the anti- 
Klein tunneling behavior is recovered. In view of the transi¬ 
tion shown and discussed above, the strength of BLG lies in 
the following fact: the tunability of bilayer graphene allows 
not only to conveniently open (and control) a band gap in its 


^ Note that an atomically sharp np junction is expected to induce intervalley 
scattering, which is not our main focus here. Numerically, such scattering 
leads to imperfection of 7 at 0 = 0 for the case of SLG, and for the case of 
BLG with ^Berry ^ as Can be noted in Fig. 9(a); see the curve marked 

by^. 



(b) 


FIG. 10. (a) Schematic of the device: a BLG flake is sandwiched 
between two h-BN layers. The whole flake can be tuned by the sili¬ 
con backgate and the central area can additionally be independently 
tuned using the central top gate. Source (S) and drain (D) are Ohmic 
contacts enabling the measurement of the conductance G. The black 
plane at the top is a schematic of the three regions deflned by the 
gates: L and R are the two leads and C is the dual-gated region, (b) 
Optical microscope image of two devices fabricated out of the same 
BLG flake. The dark yellow areas are Ohmic contacts (buried under 
BNtop) and the bright yellow ones are the top gates. The measure¬ 
ments were carried out on the right device, using the inner Ohmic 
contacts. Adapted from [24]. 
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A. Gate-tunable BLG interferometer 

The device under investigation is sketched in Fig. 10(a) and 
(b) (right device). It consists of a h-BN/BLG/h-BN stack 
which is deposited on a Si/Si02-substrate, prepared as de¬ 
scribed in Ref. [45]. The stack was realized using the dry 
transfer technique pioneered in Ref. [46] . 

For such a geometry, the silicon backgate allows us to tune 
the whole BLG stripe, whereas the local top gate acts only 
on the central region underneath. The device therefore con¬ 
sists of three areas in series, as shown in Fig. 10(a): the two 
outer regions (labeled L and R) are simultaneously tuned by 
the backgate voltage Vbg, and the central area (labeled C) 
is under the influence of both Vbg and the topgate voltage 
Vtg- The geometry is the following: the width of the whole 
flake is IT = 1.3 jim, and the lengths of the three regions are 
= 0.95 jim and ic = LI /im. 

In the following, the sample is placed in a variable temper¬ 
ature inset at temperature T = 1.6 K. We apply a constant 
symmetric bias voltage between the Ohmic contacts [labeled 
S and D in Fig. 10(a)] and record the current, allowing to mea¬ 
sure the conductance G. Additionally, the top gate voltage is 
modulated with a small AC voltage, which enables the mea¬ 
surement of the normalized transconductance JG/JVtg- 


B. Basic transport characterization 

The measurement of such a dual-gated device is shown in 
Fig. 11(a), where the conductance is displayed as a function 
of both Vtg and Vbg- Depending on the applied voltages, the 
polarity of the outer regions and of the central one can be 
changed, from hole-like to electron-like transport. This gives 
rise to four different polarity combinations: two of them ex¬ 
hibit the same polarity in the three regions (pp ’p and nn ’n - 
the prime referring to the central region C) and the other two 
exhibit opposite polarity {np’n and pn'p). The charge neutral¬ 
ity of the two outer regions is apparent from the two horizontal 
blue lines in the middle of Fig. 11(a) and the charge neutral¬ 
ity in the dual-gated region occurs along the diagonal blue 
line. The latter spans along the so-called displacement held D 
axis, as indicated in the flgure. While increasing \D\, towards 
one corner of the map or the other, the conductance at the 
charge neutrality point decreases towards zero. This insulat¬ 
ing state is due to the band gap being opened while the asym¬ 
metry between the top- and bottom layer of the BLG flake is 
made larger. Along the D-axis, the Fermi energy within the 
dual-gated region lies in the middle of the gap and the den¬ 
sity is equal to zero. As reported in various experiments [47- 
52], hopping processes resulting from residual disorder were 
found to be the dominant transport mechanism in this regime 
and at such low temperatures. 

In order to capture the electrostatic picture of the device and 
to reproduce our observations, we implemented the following 
model. First, the density n in each of the three areas - L, C and 
R - was related to the applied gate voltages Vtg and Vbg via 
the use of a parallel plate capacitor model. Each area density is 
separated into three terms: the bottom graphene sheet density. 



FIG. 11. (a) Conductance versus top and backgate voltages measured 
at 1.6 K. The purple box in the lower right corner indicates the range 
of the measurement shown in Fig. 12. (b) Calculated normalized 
conductance versus top and backgate voltages. Adapted from [24]. 


the top graphene sheet density and the intrinsic doping of the 
area, which is directly estimated from the measured data. Fol¬ 
lowing McCann [18], the asymmetry parameter u is calculated 
in each region. With the knowledge of u and n, the band offset 
is Anally calculated and inserted in the nearest-neighbor tight- 
binding Hamiltonian of BLG [9, 53]. The details of the calcu¬ 
lations are explained in Ref. [24]. As a last step, the resulting 
energy profile is used to calculate the conductance through 
the device as a function of the applied gate voltages. This 
uses a Green’s function formalism similar to Refs. [25, 28]. 
Here, only the C region is considered as the scattering region, 
L and R being treated as semi-infinite leads. As shown in 
Fig. 11(b), which displays the normalized conductance g (no 
mode counting implemented so that its maximum is 2e^/h due 
to the valley degeneracy of the spin-independent calculation), 
the overall electrostatic picture is very well captured by the 
model. 


C. Interference pattern 

As shown in Fig. 12(a), an oscillating conductance is ob¬ 
served in the pn ’p bipolar regime. To increase the visibility of 
the oscillations, one looks at the corresponding transconduc¬ 
tance map, recorded simultaneously and shown in Fig. 12(b). 
Cuts within these two maps are shown in Fig. 12(c)-(d). They 
allow to better see the strength of the oscillating pattern. In 
the following, we focus on the transconductance signal to an¬ 
alyze the oscillatory pattern in more detail. However, the same 
study could be carried out utilizing the conductance. 

As highlighted in Fig. 12(b), the oscillations evolve par¬ 
allel to the D-axis (the slope of the D-axis is displayed as a 
black line labeled /)||). This indicates that they arise from a 
mechanism taking place in the dual-gated part of the device. 
To confirm the ballistic origin of the observed signal, the fre¬ 
quency of the oscillations was analyzed. To do so, each top 
gate voltage point was converted into a density and then into 
a wavevector. Performing a discrete Fourier transform, we 
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FIG. 12. (a) Conductance measured as a function of top- and back- 
gate voltages in the pn'p regime [this region is highlighted with a 
purple rectangle in Fig. 11(a)]. Conductance oscillations are clearly 
visible, (b) Normalized transconductance map, measured in the same 
range as in (a): the oscillations appear more clearly. (c)/(d) Cuts 
along the top gate voltage axis taken from (a)/(b) (black lines) at 
Vbg = “26 V. Adapted from [24]. 


found that the oscillation frequency was A = 2.2 /im. This 
confirmed that the main contribution to the phase arises from 
interference between a directly transmitted wave and a wave 
which is transmitted through the cavity after bouncing once 
back and forth in the cavity. This corresponds to a phase dif¬ 
ference A<l> = k • 2lc, with a frequency A = 2lc = 2.2 fim. 
With this frequency analysis, we convince ourselves that the 
observation is related to ballistic transport in the dual-gated 
region, leading to a clear FP interference pattern. 


D. Broken or not broken? 

From the previous analysis, a question remains to be ex¬ 
plored: does the band gap play a role in our experimental ob¬ 
servation? To elucidate this question, we first compare the 
measured oscillations shown in Fig. 12(b) with the theoreti¬ 
cal predictions for two different systems: on one hand, a BLG 
system where the asymmetry between the layer is ignored is 
considered (i.e. the process of band gap opening is ignored 
and the gates are only responsible for tuning the overall sheet 
density) and, on the other hand, the same calculation method 
as already implemented in Fig. 11(b) is used, which follows 
McCann’s model [18] and takes into account the asymmetry. 






(b) 


FIG. 13. (a) Calculated normalized conductance as a function of top- 
and backgate voltages in the case for gapless BLG (left) and gapped 
BLG (right) (partly adapted from [24]). The latter considers the gate- 
tunable asymmetry u following Ref. [18]. (b) A line trace cut from 
the maps of (a) indicated by the dashed lines therein. 


The results of both calculations are shown in Fig. 13(a)-(b). 
Qualitatively, we observe that the shape of the oscillations is 
very different: in Fig. 13(a) (left panel), the oscillations ap¬ 
pear very sharp and asymmetric, which is not the case in right 
panel, where the oscillations appear more symmetric, almost 
like a sine function [also visible in Fig. 13(b)]. Similar obser¬ 
vations were made in Ref. [44] . Comparing these maps with 
the measurement data shown in Fig. 12(b), we conclude that 
our observation is qualitatively closer to the case where the 
band gap is implemented. 


Another confirmation lies in the visibility of the interfer¬ 
ence pattern. The visibility of the measured conductance os¬ 
cillations is found to be V = AG/Gmean = 1.5%. In a fully 
ballistic situation as those shown in Fig. 13(a)-(b), the visibil¬ 
ity of the calculated conductance yields v = 38% for the case 
where u is implemented and v = 1.5% when u is set to zero. 
However, the assumption of a fully ballistic device has to be 
weakened. Indeed, by comparing the maps shown in Fig. 12 
and Fig. 13, one can see that the measurement signal exhibits 
imperfections. We would therefore expect the theoretical visi¬ 
bility values to be upper bounds compared to the experimental 
signal. We therefore conclude that the value v = 1.5% pro¬ 
vided by the non-gapped case is inappropriate to describe our 
observation, pointing towards a role played by the band gap 
on our ballistic interference signatures. 
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FIG. 14. (a) Calculated Berry phase as a function of top- and back- 
gate voltages within the dual-gated area C of the device, (b) Zoom-in 
of the bipolar block indicated by the white box in (a). Adapted from 
[24]. 


in the map), we cannot probe the dramatic change of the Berry 
phase value. The values taken by the Berry phase in the re¬ 
gion of interest are shown in Fig. 14(b). We see that the Berry 
phase evolves, from 1.22;r to 1.46;r and is not constant. The 
limited evolution leads to only very small offsets of the os¬ 
cillations as a function of top gate voltage, such that it is not 
visible in the experimental data. Further investigations would 
be required to probe the effect of the Berry phase, ideally with 
oscillations visible in a broader range of gate voltages. 

Moreover, turning a BLG pn'p into a SLG-like pn'p by ma¬ 
nipulating the gap and hence the Berry phase is theoretically 
possible, as shown in Section III. This would however require 
a system with three cavities, each being dual-gated. Addi¬ 
tionally, the probing regime would have to cover the Berry 
phase of n (which is supposed to be close to the gap). This 
would constitute one big step beyond the above-presented ex¬ 
periment where only one dual-gated region was designed. 


E. Gate-tunable Berry phase 

As explained in Section III, signatures of the broken chi¬ 
rality due to band gap opening should appear in the Berry 
phase, which varies as a function of the induced asymmetry 
and might be accessible in experiments. To find this out, the 
oscillations have to be investigated at varying magnetic field. 

Based on the phase difference between a transmitted and a 
twice reflected electron wave, the FP resonance condition is 

^WKB + ^AB + ^Berry = j ^ (10) 

As explained earlier in Sec. IIC, <I>wkb and <I>ab are 
magnetic-field-dependent and are the origin of the parabolic 
trend of the oscillations evolving as a function of B. What 
strongly differs from the previously mentioned case of SLG 
is the Berry phase. In SLG, due to perfect transmission at 
normal incidence, the system requires finite magnetic field to 
build up trajectories which, in momentum space, enclose the 
origin and therefore pick up the Berry phase. In the case of 
gapped BLG, since even at 5 = 0 the normal incident trans¬ 
mission is finite (in the studied energy range with finite u), 
there already exist trajectories which go through k = 0, im¬ 
plying that the Berry phase is already involved.^ The Berry 
phase therefore does not depend on the magnetic field [24], 
but only on the asymmetry u. 

Figure 14(a) shows the predicted Berry phase of the dual¬ 
gated region of the device in the range of voltages available. 
We see that the Berry phase varies from 0 to 271 in this map, 
and especially involves the value 71 (white area), which is 
characteristic of SLG. However, since the interference is only 
visible in a limited range of gate voltage (white square marked 


^ Recall the trajectories sketched in Fig. 5, which do not exist in the case 
of SLG due to the Klein tunneling. Here for gapped BLG, due to finite 
transmission and refiection at normal incidence, such trajectories do exist, 
so that the ^ = 0 origin is always enclosed, independent of B. 


V. CONCLUSION AND OUTLOOK: EXPERIMENTS 
BASED ON THE BROKEN CHIRALITY 

We have explored theoretically the chirality of charge car¬ 
riers both in SLG and BLG. We pointed out experimental re¬ 
alizations of FP interferometers which allow, in the SLG case, 
to directly probe the key parameter which is the Berry phase. 
We next focused on single- and bilayer graphene systems in 
which the inversion symmetry is lifted. Such systems exhibit 
band gaps in their band structure. Close to the edge of the gap, 
the pseudospin is found to be fully z-polarized, indicating a 
broken chirality, which is then restored going away from the 
gap. This is associated with a Berry phase varying from zero 
to its pristine value (;r for SLG or 271 for BLG). Finally, we 
focused on the breaking of the chirality in BLG. Beyond the 
tunability of the anti-Klein tunneling, one of the key results 
is that the physics of BLG can mimic the behavior of SLG as 
the tunability of the transmission function allows to recover a 
perfect transmission at normal incidence for a certain energy 
range, together with a Berry phase of 7t, characteristic of the 
Klein tunneling in single-layer graphene. We presented an ex¬ 
periment on a dual-gated BLG graphene device where Fabry- 
Perot interference were probed, indicating ballistic transport. 
This allowed to illustrate the chirality breaking in more detail. 

Very recently, experimentalists have been able to use the 
pseudospin degree of freedom both in SLG and in BLG. In 
2014, Gorbachev et al. [37] demonstrated the generation of 
topological currents in a SLG system with broken inversion 
symmetry. To do so, they aligned the graphene flake on a 
hexagonal boron nitride substrate, inducing an imbalance be¬ 
tween the A and B atomic sites. Because graphene has two 
valleys, a broken inversion symmetry results in the creation of 
topological currents with different signs in each valley. This 
was confirmed by the observation of a non-local signal, as 
strong as the applied current and at micron distances from its 
path. The same was observed in a similar geometry in dual¬ 
gated BLG [51, 52]. This time, dual-gating was responsible 
for breaking the inversion symmetry. A non-local resistance 
was measured, surviving up to high temperatures. Such ad- 
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